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Abstract: Phylogenetic relationships of 11 bumblebee species, including 5 subgenera: Bombus (5 species), Tho- 
racobombus (3 species), Mendacibombus (1 species), Fervidobombus (1 species) and Pyrobombus (1 species), were 
analyzed based on the 357 bp mitochondrial cytochrome b gene sequences. There are 65 singleton polymorphic sites and 
71 parsimony informative polymorphic sites in this DNA segment, and 45 polymorphic sites within the total 119 translated 
amino acids segment. Both NJ tree and MP tree show that Mendacibombus ( B . avinovielllus ) is basal to others, followed 
by Fervidobombus ( B . pensylvanicus); Pyrobombus ( B . impatiens) and Bombus are sister subgenera; the subgenus of 
Bombus is monophyletic, in which B . ignitus diverged first. 
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(1. thESEEAE Bg. EME ABLE BE HELA RUE SORT ВМА ЗСД, ЕҢ 200025; 
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WE: Бі 357 bp HARER b НУЛЕ], Т ГЕЈ 5 11 RRR RR KAR. 5 ERE 
ТЕ Bombus (5 FR). Thoracobombus (3 FR). Mendacibombus (1 fh). Fervidobombus (1 RR) FA Pyrobombus (1 
A. BRIA 65 T AE SE Su АЛП 71 PHAR ESA, BIPERE 119 SRR A 45 THERE 
SÉQA. RU P- BEAM BEM (NJ tree) MRAMAM (МР tree) Aiba |8) А45 9: Mendacibombus 
(B . avinovielllus) 4115 8 А; Fervidobombus (B . pensylvanicus) WÈ; Pyrobombus ( В. impatiens ) Ж Bombus № 
НАШ ЕЕ; Bombus AVE RES ЖШ. ЖР В. ignitus ERWAN 5 T Ripe. 


XS: BESER AEA b EA; DNAS; GENIS; TARDES 
Ф552 =: 0969.54; 0951.3; 075 XRRR: А MBAS: 0254— 5853(2002)05 - 0361 - 06 


There are 239 bumblebee species recognized in (Williams, 1991). So some researchers have turned 
the world ( Williams, 1998), but the phylogenetic re- to molecular markers (Stephen & Cheldelin, 1973; 
lationships among them are still poorly understood be- Pekkarinen et al., 1979; Pamilo et al., 1987; 
cause of their “monotonous” morphology ( Michener, Pedersen, 1996; Koulianos, 1999; Koulianos & 


1990), especially within subgenus, such as Bombus Schmid-Hempel, 2000). However, in Bombus s. 
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1987; Ped- 
ersen, 1996) did not concur with morphology results 
(Plowright & Stephen, 1973; Williams, 1994). The 
phylogenetie relationships of those bee species are dif- 
1976; 
& 


str the molecular results (Pamilo et al., 


ferent from author to author 
Williams, 1994; 
Schmid-Hempel, 2000). 

Mitochondrial DNA (mtDNA ) is a good genetic 


marker because of its high mutation rate and haploid 


( Sakagami, 
1996; 


Pedersen, Koulianos 


without recombination (Brown et al., 1979; Avise et 
al., 1987; Moritz et al., 1987). The honeybee’s 
complete mtDNA genome was sequenced (Crozier & 
Crozier, 1993), which makes the genetic analysis of 
bees much easier. 

There are 10 species in subgenus of Bombus , and 
1998). 


Here we sequenced cytochrome b gene to study the 


most of them distribute in Asia ( Williams, 


phylogenetic relationships among 11 species including 
subgenus Bombus and four other subgenera with Apis 


mellifera as the outgroup. 


1 Materials and Methods 


1.1 Materials and PCR sequencing 

Bumblebees were collected from Beijing ( Bombus 
ignitus and B . hypocrita ) , Jilin ( B. patagiatus) and 
Zhejiang ( B . ignitus) in China, Korea ( B . ignitus ) 
and Japan ( B . ignitus ) during 2000 and 2001. The 
samples were kept in pure alcohol at room tempera- 
ture. Total DNA from 16 B .ignitus, 4 B .hvpocrita 
and 2 B.patagiatus was extracted from the heads, 
thorax (Gamery et al., 1991) or legs ( Koulianos, 
1999) using the standard phenol/chloroform methods . 

Partial cytochrome b gene amplified with the 
primer СВІ апа CB2 (Simon et al., 1994) coordi- 
nates 11 400 — 11 425 and 11 859 - 11 884 as de- 
scribed in Crozier & Crozier (1993). The reaction 
mixture consisted of 10 — 20 ng total DNA, 2 uL mixed 
dNTP (4 mmol/L), 1 pL each kind of primers (10 
pmol/uL), 3 pL MgCl (25 mmol/pL), 5 pL buffer 
(10 x, 500 mmol/L KCl, 100 mmol/L Tris-HCl, 
pH 9.0, Triton X – 100), 2 uL BSA (bovine serum 
albumin) and 2 U Taq DNA polymerase, then adding 
ddH,0 to total 50 uL. After an initial denaturation at 


96 *C for 3 min, 35 cycles were performed with dena- 
turing at 96 “С for 1 min, annealing at 54 % for 1 min 
and extension at 72 C for 1 min, followed by an end 
extension for 5 min at 72 °С. The PCR products were 
sequenced with BigDye™ Terminator cycle sequencing 
kits ( Perkin-Elmer Biosystems ) using ABI377 automat- 
ie sequencer (PE). Each individual was sequenced in 
both directions. The sequences were checked by eyes 
and edited with DNAstar (DNASTAR Inc.1996) pack- 
age. Cytochrome b gene sequences GenBank acces- 
sion numbers of B.hypocrita, B . ignitus and B . pat- 
agiatus are AF 498397, AF498398 and AF498399. 
1.2 Sequences from GenBank 

The sequences in 8 bumblebee species of 5 sub- 
genera and Apis mellifera were collected from GenBank 
(Table 1). 


Table 1 Sequence data from Genbank in bumblebee 
and honeybee species 








Subgenus Species No. of sequence 

Bombus B. lucorum AF002722 

B . terrestris AF002721 
Thoracobombus B . pascuorus AF017516 

B . humilis AF017517 

B . ruderarius AF002723 
Mendacibombus В . avinoviellus AF181610 
Fervidobombus B . pennsylvanicus AF181611 
Pyrobombus B . impatiens AF281169 

Apis mellifera L06178 


1.3 Phylogenetic analysis 
The polymorphic sites and synonymous/non-syn- 
substitutes software 


onymous were calculated by 


DNASP 3.0 (Rozas & Rozas, 1999). We use the 
amino acid sequence to construct the phylogenic trees. 
NJ tree was constructed with MEGA version 2.1 (Ku- 
mar et al., 2001) based on P-distance (Nei & Ku- 
mar, 2000), while maximum parsimony tree with 
PAUP*4.08b (Swofford, 2001) by heuristic search 
with TBR branch swapping. Apis mellifera was used as 
the outgroup. For both NJ and MP tree, we used 
2000 bootstrap replications to estimate the support 
rate. We reconstructed NJ tree of subgenus Bombus 
based on DNA sequence to further analyze the phyloge- 
ny, because within subgenus DNA sequence can re- 


flect the relationship better than amino acid sequence. 
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2 Results and Discussion 


There are 65 singleton polymorphic sites and 71 
parsimony informative polymorphic sites in this DNA 
segment; 45 polymorphic sites within the total 119 
translated amino acids segment; but no gap or terminal 
code were found (Fig.la, b). The mutation pattern 
of the segment suggests that it is probably from mito- 


chondrial genome, not from nuclear psuedogene. Be- 


cause the mean transition/transversion is 0.6347 and 
the rate of amino acid polymorphic site is high 
(37.82%), we used the amino acid sequence to con- 
struct the phylogenetic trees. Both NJ and MP trees 
showed the same pattern: Mendacibombus (В . avino- 
vielllus ) is basal to others, followed by Fervidobombus 
(B . pensylvanicus ); Pyrobombus ( B . impatiens ) and 
Bombus are sister subgenera (Fig.2a, b). The sub- 


genus of Bombus is monophyletic: B . locrum , B . pat - 
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Fig.1 Segment mutation sites of cytochrome b DNA and amino acid sequence in 12 bumblebee 


and honeybee species 
a. 357-segment cytochrome b DNA sequence; b. 
1) B. impatiens, 2) B.avinoviellus, 3) В. pascuorus 


119-segment cytochrome b amino acid sequence. 
, 4) B.humilis, 5) B.locrum, 6) В. pennsylvanicus, 7) 
B . terristris, 8) B.ruderarius, 9) B . ignitus, 10) B. 


hypocrita, 11) B.patagiatus, 12) Apis mellifera . 


Figure above those sequence sites is the site position of the segment sequence; dotes indicate the nucleotides identical to 


that of В. impatiens . 


agiatus and В .hypocrita are closely related to each 
other, and B.ignitus diverged first among those five 
species (Fig.2a, b). MP tree length is 83, Consis- 
tency Index is 0.819, Retention Index is 0.712 and 
Rescaled Consisteney Index is 0.583. Similar results 


for the subgenus of Bombus were obtained using DNA 


sequence directly, except that B.patagiatus and 


B . hypcrita become the most derived sister species. 
The phylogeny of bumblebees has not been re- 

solved yet because of “ monotonous” morphology 

(Williams, 1991). Now by the analysis of the mito- 


chondrial eytochrome b sequences , we provided new 
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B.ignitus 
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B.humilis | Thoracobombus 
68L B ruderarius 
| Fervidobombus 
| 
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B.impatiens Pyrobombus 


B.lucorum 
77 


B.hypocrita 


61 Bombus 


B.patagiatus 
B.terrestris 
B. ignitus 


B.pascuorus 


93 


B. humilis Thoracobombus 


B.ruderarius 





B.pennsylvanicus | Fervidobombus 





B.avinoviellus і Mendacibombus 


Fig.2 Phylogenetic trees of 11 bumblebee species based on mitochondrial cytochrome 
Ь gene sequences 
a. NJ tree constructed with MEGA V2.1 based on P-distance; b. Topology from MP analysis using 
PAUP” V4.08b. 
The numbers at the nodes are bootstrapping indices of 2 000 replicates of the level of support for individual 
nodes, and data are dismissed when the support rate is less 50%. The corresponding subgenera are 


shown at the right. 


information on phylogenetic relationship of bumble- and Pyrobombus were paraphyletic between two 
bees. clades, but the latter ones showed that Bombus was 

The topology of subgenera from our molecular data basal to Pyrobombus in a monophyletic clade. The re- 
is similar to that of Williams! morphologic study sult based on cytochrome b sequences was similar to 


(Williams, 1994), but differs from other recent that of Pedersen (1996) for the two subgenera ( Kou- 
molecular studies. Pedersen ( 1996) showed that lianos & Schmid-Hempel, 2000). However, all the 


Bombus was basal to Pyrobombus and several other studies shared some similar patterns that the subgenus 
subgenera on the basis of COI sequence. Even based Bombus is monophyletic . 

on the same COI sequences, the results of Koulianos Because of the high rate of polymorphic sites a- 
& Schmid-Hempel (2000) differed from that of Peder- mong our analysis data, we believe that mutation of 
sen (1996). The former results indicated that Bombus mtDNA among subgenera is almost saturated, so the 
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translated amino acid sequence is better than DNA se- 
quence for the phylogenetic analysis of bumblebee sub- 
genera. In the future, it is necessary to study more 


species using other more conservative genes. 
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